Lower bound for electron spin entanglement from beamsplitter current correlations 
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We determine a lower bound for the entanglement of pairs of electron spins injected into a meso- 
scopic conductor. The bound can be expressed in terms of experimentally accessible quantities, the 
zero-frequency current correlators (shot noise power or cross-correlators) after transmission through 
an electronic beam splitter. The effect of spin relaxation (Ti processes) and decoherence (T2 pro- 
cesses) during the ballistic coherent transmission of the carriers in the wires is taken into account 
within Bloch theory. The presence of a variable inhomogeneous magnetic field allows the determi- 
nation of a useful lower bound for the entanglement of arbitrary entangled states. The decrease in 
entanglement due to thermally mixed states is studied. Both the entanglement of the output of a 
source (entangler) and the relaxation (Ti) and decoherence (T2) times can be determined. 



Quantum nonlocality has been an intriguing issue since 
the early days of quantum mechanics Q . Nonlocal effects 
can come into play when a quantum system is composed 
of at least two subsystems which are spatially separated. 
Despite their simplicity, the Bell states of two distant 
quantum two-state subsystems (A and B) 
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exhibit the essential phenomenology of quantum nonlo- 
cality (e.g., they violate Bell's inequalities 0) thus pro- 
viding an ideal testing ground for quantum nonlocality. 
Here, we represent the two-state systems as spins 1/2 
with basis states "spin up" 1 |) and "spin down" 1 1) with 
respect to an arbitrary fixed direction in space. 

With the development of quantum information theory 
and in particular with quantum communication, it 
has become clear that EPR pairs can also play the role of 
a resource for operations that are impossible with purely 
classical means. In this context, two-state systems are 
referred to as quantum bits (qubits), and quantum non- 
locality is related to the concept of entanglement (defined 
below). A number of quantum information processes- 
quantum teleportation [3 , quantum key distribution [j| , 
quantum dense coding |f|, etc. -have been successfully im- 
plemented using pairs of photons with entangled polar- 
izations, i.e., in states such as Eqs. Q and (J2J). Pho- 
tons have the advantage that they are easily moved form 
one place to another, allowing for experiments involving 
space-like separations between detection events 0. 

More recently, there has been increasing interest in the 
use of the spin of electrons in a solid-state environment 
for spin-based electronics Q and as qubits for quantum 
computing ||. Subsequently, quantum communication 
on a mesoscopic scale, typically on the order of microm- 
eters in semiconductor structures (e.g. quantum wires), 
was proposed ||. Rather than achieving space- like sep- 
aration between detection events on the two sites (this 




FIG. 1: Inset: Proposed setup with two-electron scattering 
at a beamsplitter (BS) with transmittivity T. Electrons are 
injected pairwise from the entangler (E) into contacts 1 and 
2 of the BS. The mean current / = (I a ) and one of the cur- 
rent correlators S a p are measured at the outgoing contacts 
a,P — 3, 4. Plot: Entanglement of formation E of the elec- 
tron spins versus singlet fidelity F and the reduced correlator 
/ = S 3S /2eIT(l - T). The curve illustrates the relation be- 
tween noise and entanglement for Werner states. For general 
states, the curve represents a lower bound for the entangle- 
ment, i.e. allowed values for E and / (or, equivalently, F) 
are represented by points in the shaded region. Measuring / 
determines a lower bound for the entanglement E. 



would require sub-picosecond detection), the idea here 
is to use quantum entanglement between parts of a co- 
herently operating solid-state device (in the most ex- 
treme case, a quantum computer). It is then relevant 
to study the transport of spin-entangled electrons in a 
many-electron system and possible means of entangle- 
ment detection. Two-particle interference at a beam- 
splitter (BS) combined with the measurement of current 
fluctuations (in general, the full counting statistics 
[l"lj|) was identified as a detector for entanglement. 

In this paper, we go one step further, providing a lower 
bound for the amount E of spin entanglement carried by 
individual pairs of electrons, related to the zero- frequency 
current correlators when measured in a BS setup (Fig.^ 
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Inset) by injecting the electrons separately into the two 
ingoing leads (1 and 2) and measuring either the current 
autocorrelator S aa in one of the outgoing leads (a = 3, 4) 
or the cross-correlator S34. It is assumed that the size of 
the scattering region is smaller than both the coherence 
length and the mean free path, allowing for ballistic and 
coherent transport. In the following, T will denote the 
transmittivity of the BS, i.e. the probability to be scat- 
tered from lead 1 to lead 4 (or from 2 to 3). The ideal 
BS for the proposed setup does not give rise to backscat- 
tering (e.g. from lead 1 back into lead 1, or from 1 into 2, 
etc.). We will also analyze the effect of such backscatter- 
ing processes, as they give rise to background shot noise 
which is unrelated to entanglement. During their trans- 
port, the electron spins will be exposed to decoherence 
and relaxation due to spin-dependent scattering caused 
by magnetic impurities, nuclear spins, or the spin-orbit 
coupling (see [12( for a review). We include these effects 
within a Bloch equation formalism |l8l | . Comparison be- 
tween our theory and experiment will (i) test proposed 
entanglers Q EE E E3 El G3 El and (ii) determine 
spin relaxation (T±) and decoherence (T2) times. 

The materials and structures required for testing our 
theory, although at the forefront of current capabilities, 
appear to be feasible. The largest efforts seem to be 
necessary to realize the electron spin entangler |lOj for 
which there exists a number of theoretical ideas, us- 
ing norm al- 1 14l Il5ll o r carbon-nanotube-superconductor 
junction s Il6lll7l ll8 | . or single E3|j ° r coupled quantum 
dots (iTtL l20l| . The electronic BS and the measurement of 
BS current correlators have been experimentally demon- 
strated in a GaAs/AlGaAs heterostructure |2lJ. Coher- 
ent transport of electron spins over more than 100 pm in 
GaAs has been observed |22|. 

Traditionally, current correlations, and in particular 
the quantum partition (shot) noise have been used to 
gain information about a scatterer beyond its conduc- 
tance |23|. Here, we use a known scatterer (the BS) to 
gain information about the quantum state (more pre- 
cisely, its entanglement) of the scattered particles. The 
correlation function between the currents I a (t) and Ip (t) 
in two leads a, (3 — 1, ..,4 of the BS is defined as 



hv f T 
S a n(u) = lim — / dt 
7-^00 r J 



ReTY[SI a (t)6I (O)x\, (3) 



where SI a = I a — (I a ), (I a ) — Tr(/ Q x); v 1S the density 
of states in the leads, and X is the density matrix of the 
injected electron pair (below, we suppress the orbital part 
of X, see [HJ for Coulomb effects). Writing x in the Bell 
basis, Eqs. Q and ©, x = +G |* + )(*+| + 

J2i=± G» I I) and S a p = S a p(uj = 0), we arrive at 

S a , = FS^ + G S^ + £ Gi S%\ (4) 

»=± 

S^J ee lim — [ T dtRe(y\6I a (t)5I (O)\*). (5) 



Using the standard scattering approach [23j, we have 
found earlier E3 that the singlet state |^/_) gives rise 
to enhanced shot noise (and cross-correlators) at zero 
temperature, 4*^ = = 2eiT(l - T)f, with 

the reduced correlator f = 2, as compared to the "clas- 
sical" Poissonian value / = 1 |24| . The average currents 
are given by I = (I 3 ) = (I4) = e/hv. We also know 



that all triplet states are noiseless, S a 







S. 
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(a, (3 = 3,4). Both the current autocorrelations (shot 
noise) and cross-correlations are only due to the singlet 
component of the incident two-particle wavefunction, 

S33 = —S34 = FS\ 9 ~) = 2eiT(l — T)f, f = 2F. (6) 

Including backscattering with probability Rb, we find 

5 33 = 2eI[2F(l-R B )T(l-T)+R B /2], (7) 

5 34 = -2e/2F(l-i? B )T(l-T), (8) 

where / = (e/hv){l-R B ). Since/' = S 3i /2eIT(l-T) = 
2F(1 — Rb) < / is smaller than / without backscatter- 
ing and the entanglement of formation E is a monotonic 
function of / (see below and Fig. QJ, we still obtain a 
lower bound on E (the bound will become less informa- 
tive as Rb increases). Note that this does not hold for the 
autocorrelator S33. However, one can determine Rb, e.g. 
by measuring the shot noise power using normal Fermi 
lead inputs [2l| and then obtain / from either 633 or S34. 

The entanglement of a bipartite state x G 7~La ® 
can be quantified by its entanglement of formation |2J| 
E(x) = min{(|xi),ft)}e£(x) Y, t Pt s N(\Xi}), where £{x) de- 
notes the set of ensembles {(|Xi)>Pi)} f° r which x = 
^2iPi\Xi)(Xi\- We have used the von Neumann en- 
tropy of the reduced density matrix ps = Trx|V')(V'l) 
Sn(\iP)) = — Tr^/Os logps (logarithms are in base 2). 
A state with E > (E = 1) is (maximally) entangled, 
whereas a state with E = is separable (in the case of 
a pure state, it is a product i/ja ® ^>b)- The Bell states 
Eq. and J5J are maximally entangled. Neither local 
operations nor classical communication (LOCC) between 
subsystems A and B can increase E. In quantum infor- 
mation theory, E(x) is the maximal ratio N/M of the 
number N of EPR pairs (maximally entangled states) 
required to form M copies of x as A — > 00; E is the 
quantity that measures how much of the resource (quan- 
tum entanglement) is available. 

For arbitrary x, E(x) cannot be expressed as a func- 
tion of only its singlet fidelity F = However, 
this is possible for the so-called Werner states [26j 



p F = F|*_)(*. 



l-F , 



l* + ><* + l+£l^X^I 



i=± 



(9) 

being the unique rotationally invariant states with sin- 
glet fidelity F. It is k nown 25] that E(F) = E(p F ) = 
H 2 (l/2 + y/F(l - F)) if 1/2 < F < 1 and E(F) = 
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FIG. 2: Homogeneous magnetic field Sh — and P = 1. 
(a) / of the spin singlet state |^-) after ballistic transmission 
through a BS as a function of the spin decoherence time T2, in 
units of the ballistic transmission time to = L/vf- Different 
curves correspond to different values of the spin relaxation 
time Ti (same units). Note that we only plotted the curves 
for T 2 < 2Ti. (b) / for a spin triplet state |*+). Since / < 1, 
the lower bound on entanglement is zero, i.e. we cannot learn 
anything about entanglement of injected triplets at Sh — 0. 
(c) Lower bound on the entanglement of formation E. 



E(p F ) = if < F < 1/2, with the dyadic Shannon 
entropy H2{x) — — xlogx — (1 — x) log(l — x). Together 
with Eq. ©, this enable us to express the entanglement 
of Pf in terms of the reduced correlator / (Fig.^l. 

We generalize this result to arbitrary mixed states \ 
of two spins (qubits). Any state x can be transformed 
into the Werner state pf with the same singlet fidelity F 
by a random bipartite rotation [251 \21\. i.e. by choosing 



a random U € SU (2) and applying U ® U to X- Since 
this operation involves only LOCC, E cannot increase, 



E{F) < E( X ). 



(10) 



The entanglement of formation E(F) of the correspond- 
ing Werner state therefore provides a lower bound on the 
entanglement of x (FigCJ. Thus, a noise signal exceeding 
the Poissonian limit (/ > 1) in the BS setup can in prin- 
ciple be interpreted as a sign of entanglement between 
the electron spins injected into leads 1 and 2 pjof . 

We now include relaxation and decoherence into our 
analysis. At time t = 0, we start with a spin singlet 
(upper sign) or triplet (lower sign) state 



*(0) = |* T )<* T |. 



(11) 



We describe the dynamics of x(t) m & field B || z and in 
the presence of spin decoherence (T 2 processes) and relax- 
ation (Ti) phenomenologically within a single-spin Bloch 
equation for the polarization P = ((a x ), (o- y ), (<t 2 )), 



P = P x h R{P - P) = -Q(P - P) 



(12) 



with (<7i) = Tr(CTip), h = gps^ = (0,0, h), the station- 
ary polarization P = (0, 0, P) (note that Pxh = 0), and 
the relaxation matrix |3l| Rij — SijRi with Ri = i? 2 = 
77 1 and R 3 = Tf 1 . Solving Eq. we obtain 



P(t) = e- s "P(0) + (1 



-m 



)P, 



or, in terms of the spin density matrix, 

p(t) = (P +P(t)- < T)/2 = A h (t)[p(0)}, 
with the superoperator (a(t) = 1 — e -t / Tl ) 



(13) 



(14) 



5(PTT + Pll)(l + "(ft?) + l(Pn ~ Pll)e- t/Tl 



Mm 



-t/T 2 +iht 



s(ptt + pu)(i - a ( f ) p ) - |(ptt - Pn) e 



-t/T! 



r 



(15) 



with the matrix elements pij = (i\p\j) and j =T,4)- 
We apply A^(i) to both spins individually, 



X (t) = (A hl (t)®A h2 (t))[ X (0)}, 



(16) 



where hi is the field at electron i. Using Eq. © and 
F(t) = (^_|x(i)|*_) at time t > 0, we obtain 

f(t) = ±e" 2t/T2 cos(Sht) 

+ i(l + e- 2 */ Tl ) - -(1 - e-'^fP 2 , (17) 

(1 2) 

where 5h = hi — h^. If the decoherence time T 2 ' of the 



two electrons is different, then T 2 in Eq. (|17f) becomes 
T 2 EPR = (1/T 2 (1) + 1/rfV 1 - We define Tf PR similarly 
if P = 0. However, if P = 1, then exp(— t/T{) is replaced 
byexp(-i/T 1 (1) ) + exp(-i/T 1 (2) ). 

A homogeneous field, Sh = 0, docs not affect /. For 
slow relaxation, T x > t, we find f(t) = 1 ± e" 2t/T2 . In 
Fig. we plot / for Sh = and P — 1 versus T 2 in 
units of the ballistic transmission time [24| to = £/Vf 
(L=length of ballistic trajectory, Wi?=Fermi velocity). 
For unentangled triplet states, x(0) = | TTXTT I, I 44X44 I, 
we find / < 1/2 for all T u T 2 , and P (Fig.Eb). 
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(Fig. 03). Above the critical temperature T c = 0.91 J/fcs 
there is no entanglement even for T\,Ti — > oo. 
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FIG. 3: (a) Reduced current correlator / versus field inho- 
mogeneity 5h = hi — hi (in units of h/tog^B, io=ballistic 
transmission time, p=g-factor, /is=Bohr magneton) for an 
injected singlet state and P = 1. Solid lines represent Ti S> 1 
and T2 = 1, 3, 10, 00, grey dashed lines Ti = 3 and T2 = 1, 3. 
For injected triplets, the plot is phase-shifted by n, providing 
tight lower bounds at Sh = tt. Tight lower bounds for any 
input state can be determined by varying the direction of Sh. 
(b) Plot of / for a thermally mixed initial state versus T2 (in 
units of to) for Ti ^> t and P = 1. The various curves corre- 
spond to k B T/J = 0, 0.25, 0.33, 0.5, 1, where J and T denote 
the exchange energy and temperature during the preparation 
of the state. Inset: The maximal / (at T\,T^ to) versus 
fcsT/J. There is no entanglement (/ < 1, E = 0) above the 
critical temperature T c = 0.91 J/ks- 



An inhomogeneous field Sh ^ (or, equivalently, a lo- 
cal controllable Rashba spin-orbit coupling [2£|) has the 
effect of continuously rotating singlets into triplets and 
vice versa (Fig. [3J. This a lower bound of the triplet 
entanglement, x(0) = |^/ + )(^/ + |, which is as tight as 



Eq. 02J for the singlet, 



E > max£(/ (Sh) /2) . 

5h 



(18) 



where f(Sh) is the measured noise power (or cross- 
correlator), E(F) is the entanglement of the Werner state 
Pf (Fig. If a field inhomogeneity 8h can be created 
pointing in arbitrary directions in space, then the above 
result represents a tight lower bound for any injected en- 
tangled state. In particular, each maximally entangled 
state I \P) will be detected in this way, since there exists a 
u = cxp(-i<5h • er) S SU(2) such that u f <g> u|*) = |*_). 
This rotation can also be done unilaterally, i.e. there is a 
v G SU(2) with v ® = |*_) (see also [H). 

Finally, we study the case where the spin state of the 
injected pair of carriers Eq. Qllj l is mixed, because it is 
prepared at a temperature T comparable to the energy 
splitting between spin states, typically (if the Zeeman 
effect is negligible) the exchange energy J, i.e. the singlet- 
triplet splitting. In this case, x(0) = pp with F = (1 + 
3e -J / feflT ) -1 where ks is Boltzmann's constant. We only 
show the resulting / for T\ ^> t here (the full expression 
will be reported elsewhere [29|). 



/(*) 



I _|_ e -2t/T 2 + £ -J/k B T ^ _ e -2t/T 2 
J/k B T 



(19) 



l + 3e- 

which is the statistical mixture of Eq. i|17[l for the sin- 
glet and triplet with the appropriate Boltzmann weights 
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If the current is carried by quasiparticles with charge e* , 

/ will be renormalized by a factor e/e*. 

Within standard weak-coupling theory, T2 = [(2Ti) _1 + 

Tn" 1 ] -1 < 2Ti (T^ 1 is the longitudinal decoherence rate). 

The superoperator A^(t) is linear and trace-preserving 

on density matrices with arbitrary trace since we have 

not imposed the trace condition Po = 1 at this point. 



